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Statistics of heat exchange between two resistors 

D.S. Golubev and J.P. Pekola 
Low Temperature Laboratory, Department of Applied Physics, 

Aalto University School of Science, P.O. Box 13500, 00076 AALTO, Finland 

We study energy flow between two resistors coupled by an arbitrary linear and lossless electric 
circuit. We show that the fluctuations of energy transferred between the resistors are determined 
by random scattering of photons on an effective barrier with frequency dependent transmission 
probability t{oj). We express the latter in terms of the circuit parameters. Our results are valid in 
both quantum and classical regimes and for non-equilibrium electron distribution functions in the 
resistors. Our theory is in good agreement with recent experiment performed in the classical regime. 


I. INTRODUCTION 

The problem of energy exchange between two resistors 
has been first analyzed by Nyquisti on the way towards 
his famous formula for the current noise of a resistor, 

Si = AkeT/R. ( 1 ) 

Here Sj is the spectral density of noise at low frequen¬ 
cies |w| <C kBT/h, ks is the Boltzmann constant, T is 
the temperature and R is the resistance. Equation (ED 
has been confirmed by Johnsoi>^ and by numerous sub¬ 
sequent experiments. For a long time afterwards trans¬ 
port of heat in electric circuits has been considered well 
understood. Recently, however, it has attracted renewed 
attention due to advances both in theory and in tech¬ 
nology. On the theoretical side, the discovery of the 
fluctuation theorem^Ti^ has triggered the interest in the 
statistics of heat transport. Statistics of effective elec¬ 
tron temperature fluctuations in small metallic grains is 
also under discussionii^. The experiments have recently 
advanced in two directions. First, quantum transport 
of heat between two resistors coupled by superconduct¬ 
ing wires and separated by up to 50 /rm distance has 
been demonstrated at sub-kelvin temperatures^id^. Sec¬ 
ond, utilizing low noise amplifiers Ciliberto et al. have 
recently measured the full statistical distribution of heat 
transferred between two resistors kept at temperatures 
88 K and 296 K respectivel^ d^d^ . They have verified 
the validity of the fluctuation theorem and worked out a 
theoretical model based on Nyquist’s formula ED- 

Motivated by these developments, in this letter we pro¬ 
pose a theory of full counting statistics of photon medi¬ 
ated heat exchange between two metallic resistors valid 
both at high and at low temperatures, where the classi¬ 
cal formula for the noise ED can no longer be used. We 
consider two resistors, i?i and i ?2 shunted by impedances 
Zi{u!) and Z 2 {u>), and coupled by a linear element (e.g. 
transmission line, capacitor, etc.) having the impedance 
Zo{uj) (see Fig. la). The impedances Zj{oj), [j = 0,1,2) 
are purely reactive and do not generate noise. The aver¬ 
age photonic heat current flowing from the resistor 1 to 
the resistor 2 reads 

poo J 

Jq = J ^ wr(w) [m (w) - 712 (w)], (2) 
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FIG. 1. Two resistors connected in a linear circuit: (a) gen¬ 
eral case - resistors are connected by an arbitrary reactive el¬ 
ement with the impedance Zo{uj) and shunted by the reactive 
impedances (cu), ^ 2 ( 0 ;); (b) classical regime T\,T 2 > oJc, (c) 
quantum regime Ti,T 2 ^A (d) realistic model, with stray 
capacitances Gi, C 2 and wire inductances Li, 1/2; (e) resistors 
directly coupled by two ideal zero resistance wires. 


where t{u}) is the effective transmission, which we will 
specify later, rijiuj) are photon distribution functions 
(here and below we put ks = ^ = 1)- Typically t{u!) 
drops at certain cutoff frequency Wc- Assuming that 
ni(a;), 712 ( 0 ;) have equilibrium Bose form with the tem¬ 
peratures Ti and T 2 , one finds that at high temperatures, 
Ti,T 2 > utc (Fig. lb), Jq r; T(0)o;c(ri - T 2 ) in agree¬ 
ment with experimental findings of RefsJ^d^. In this 
classical regime Nyquist’s formula ED may be used to 
derive the heat current. In this letter we will be mostly 
interested in the opposite, quantum, limit Ti,T 2 ^ uJc 
(Fig. Ic), which is relevant for typical low temperature 
experiments^iiS. Indeed, the cutoff frequency may be es¬ 
timated as ojc ^ TD.\n.{\/RjCj, Rj/Lj}, where Cj ^ eeol 
are stray capacitances, Lj ^ figl are inductances of the 
wires (Fig. Id), cq and /iq are vacuum permittivity and 
permeability, e is the dielectric constant, and I is the char- 
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acteristic size of the sample. For the parameters of the 
low temperature experiments^*^, namely T ^ 100 mK, 
R I kn and I ^ 10 /im, one finds T/ojc 10“^ <C 1. 
Thus the circuit is in the quantum regime. In contrast, 
for the experiments by Ciliberto et with T ^ 100 

K, i? = 10 Mil and I ^ 1 cm one finds T/wc ^ 10^° 1, 

which corresponds to strongly classical regime. 


The trace in Eq. ([5]) can be expressed as a path in¬ 
tegral over the fluctuating potentials ,V^, 
defined on the forward (F) and backward (B) branches 
of the Keldysh contour, and over the Grassman fields 
, a^*, , a^* describing electrons. Performing the 

Gaussian integral over the latter, we get 

e^= [ VV^’^VA^’^ (6) 


II. MODEL 

Our goal is to find the distribution of the energy Q 
transferred from the resistor 1 to the resistor 2 during 
the time t, which we denote as P{t,Q). It is more con¬ 
venient to work with the cumulant generating function 
(GGF), F(l, A), which depends on the counting field A 
and defined as 


where the effective action , A^’^] is the sum of 

the electronic and electromagnetic contributions. 


iS^ = -I- iSem, 

zFii = 21n[det(G-i[y^,I/^])], 


iSoju — i 


dt' 


_ ^2 _ rr2 


Hi 


Stt 


( 7 ) 

( 8 ) 

(9) 


= J dQP^^P{t,Q). (3) 

We describe the system by a Hamiltonian 

H = Ho + Hem + Hint, (4) 

where Hq — ^kd\^dka is the Hamiltonian of non¬ 
interacting electrons moving in the combined poten¬ 
tial of ion lattice and impurities, dka- is an annihila¬ 
tion operator of an electron in the eigenstate lipka) (o 
is the spin index) and Ck is the corresponding eigen- 
energy; iJem = / d^r{E^ -I- iT^)/87r is the Hamilto¬ 
nian of electro-magnetic field; E and H are the op¬ 
erators of the electric and magnetic fields respectively; 
Hint = - J2kn,a ^^kudl^dna is the interaction Hamilto¬ 
nian; and Vkn = {'<Pk\V{f')\'4’n) are the matrix elements 
of the electric potential operator between two eigenfunc¬ 
tions of the non-interacting electron Hamiltonian Hq. 
The Hamiltonian Hq describes the two resistors, the wires 
connecting them and the leads attached to them if they 
present. 

An important point is the definition of the transferred 
energy Q. Here we have in mind the detection scheme 
based on normal metal - superconductor tunnel junc¬ 
tions attached to the resistors^d^. Such a junction al¬ 
lows one to measure the effective temperature of a resis¬ 
tor or, more generally, the distribution function, f{E, r), 
of electrons in ild^. The latter can be converted into 
the total electron energy of the resistor j {j = 1,2) as 
£j = 2 d^r J dE Evj{E)f{E,r) (here in the vol¬ 
ume of the resistor j and Fj{E) is the density of states). 
Within this approach it is natural to define the trans¬ 
ferred energy as the drop in the electronic energy of the 
resistor 1 during the time t, Q = —£i{t) -I-£ 1 ( 0 ). The 
corresponding quantum expression for the GGF reads^: 




tr 




( 5 ) 


where po is the initial density matrix and Hi is the free 
electron part of the Hamiltonian of the resistor 1. 


Here we introduced the inverse Keldysh Green function 
of electrons -k (5G):„\ where 


^-1 __ 
^0,fcn ~ 


^kn 


idt - €k 0 
0 -idt + Cfc 


AG 


-1 

kn 


0 



( 10 ) 


At this stage we retain the information about occupation 
numbers of all energy levels keeping the dependence of 
the counting filed Xk on the level index k. Below we will 
only consider linear circnits free of highly resistive jnnc- 
tions or quantum dots in the Goulomb blockade regime. 
Then one can expand the action (|S]) to the second order 
in 


This expression contains the Green function of non¬ 
interacting electrons, Gq. It is defined as 


Go,fcn(G;^2) — i^kn^ 

^ 12(1 — fk) — 02l/fe 
I- fk 




-fk 

— 0l2fk + fi*2l(l — fk) 


where 6ij = 9(ti — tj) are Heaviside functions and fk = 
{d\^dkcr) are the occupation numbers of the energy levels. 
The first term in the expansion (HU does not depend on 
Xk and may be omitted. The second term, tr [2GoAG“^], 
is canceled by a similar contribution coming from posi¬ 
tively charged ion background. Thus, only the last term 
of Eq. (fTTl) matters. We transform it to the from 





( 12 ) 


Here we have introduced the potentials V~^ = (V^ + 
V^)/2 and V~ = V^ — V^, as well as dimensionless com¬ 
binations Xkn containing electronic distribution functions 
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fk and counting fields Xk- 

= /fe(l - /n) - 1) 

+ (1 - fk)fn - 1) , 

Xtn = (012 - 02 l)(/fc - fn) + fk(l - 

Xkn = 0> 

Xku = -*(1 - fn) + 1) /4 

- (1 - fk)fn + 1) /4. (13) 

Next we perform disorder averaging of the matrix el¬ 
ements in Eq. (fT^ inside the metallic parts of the 
system ignoring weak localization and utilizing the rule of 
averaging for the product of electronic wave functionsi^ 


(V’fe(r’2)V’n(r2)V'/c(n)iA*(ri)(5(Ei - ek)S{E2 - £„)) 
kn 

= (iz/tt) ReV{Ei — E 2 ,ri,r 2 ). (14) 


Here u is the density of states and 'D{E,ri,r 2 ) is the 
solution of the diffusion equation {—iE — D{r)S/'^)'D — 
(5(ri — r 2 ), where D{r) is the diffusion constant. In 
good metals with local current-field relation, j = a{r)E, 
where a{r) = 2e‘^i^o{r)D{r) is the conductivity, one can 
approximate Re'D{E,ri,r 2 ) —)> —D{r)V'^/E'^, and the 
action cm acquires the form 


iSgi = —J dt'dt" J d^ra{r) J 


dw 

27r u) 


X ^ 77“(',VH^(t',r)VH“(t",r). (15) 

a,/3=± 

Here 

- 1) - + l)(e-*^"‘^ - 1), 

= 1 - - 1) + (llc.,r + l)(e"*^-“ - 1), 

du)^r 

+ 1) + {nu,,r + l)(e-*^’'‘" + 1) 


Vu.T = 


, ( 16 ) 


and n^j^r is the effective photon distribution function, 


1 

^io,r — 

(jJ 




(17) 


It satisfies n-cj,r = — 1 —Hw,?- and in local equilibrium, i.e. 
for momentum isotropic electron distribution function of 
the form f{E,r) = l/(e^/^*^'^^ -I- 1), where T{r) is the 
local electron temperature, it reduces to Bose function 
1/(6“^/^*^— 1). However, ,, may deviate from simple 
Bose form if the electron distribution function is driven 
out of equilibrium by, for example, bias voltage applied 
to a resistor—. In Eq. (ESI we have also assumed that 
the counting field is the same for all energy levels with 
wave functions localized in the vicinity of the point r and 
that it slowly varies in space at distances exceeding the 
spatial extension of these wave functions. 


We are now in position to write down the action of two 
coupled resistors depicted in Eig. la. We put A(r) = A^ , 
cj{r) = cjj {j = 1, 2) inside each resistor. Considering low 
frequency modes, we also put VH(r) = Vj/Lj, where Vj 
is the instantaneous voltage drop across the j—th resis¬ 
tor, and Lj is its length. We also define the resistances 
Rj = Lj/ajAj, where Aj are the cross-sectional areas of 
the resistors. With these approximations we get 


doj 


27r 


ujRj 


= -Y 

1=1.2 “'o 

X [r,YA)V+A)V+{t") + vt~A)V-{t')V+it") 

+ rj--A)V-it')V-in], ( 18 ) 


where the functions 77 “^ (w) are given by Eqs. m with 
photon distribution functions averaged over the volume 
of the resistors, njA) = Jq. d^rriuj^r/^j, and with Xr 
replaced by Xj. The fields E and H in 3d space around 
the resistors and other circuit elements can be expressed 
via the voltages Vj by solving linear Maxwell equations 
with proper boundary conditions. In this way one finds 


EF,B{t, r) = 




f dt'[ei{t-t',r)vA^{t') 

J —00 

+ e2{t-t' ,r)V2'^{t')], (19) 

f dt'[hi{t — t'A) 

J —00 

+ h2{t-t\r)V,^’^{t')], (20) 


where ej{t, r) and hj{t, r) are the fundamental solutions 
for electric and magnetic fields, which depend on the 
sample geometry. The solutions (jl9l20l) should be sub¬ 
stituted into the electro-magnetic part of the action 
After the integration over coordinates, this action be¬ 
comes quadratic in the potentials Vj. Moreover, since 
El - El - Hi + H% = 2E-E+ - 2H-H+ only the 
combinations VA'^l' appear in it. The coefficients in 
front of these combinations are expressed in terms of the 
functions ej(t, r), hj{t, r) and determine the impedances 
ZjA)-: shown in Eig. la, for a given sample. Einally the 
electro-magnetic part of the action acquires the form 




/■* 

f duj 



- / dt'dt" 




Jo 

J 27 r 

UJ 


V 

1 

+ 

in , 

V2(l')Vt(l")l 


'- 1 = 1.2 

ZjA) 


ZqA) 


( 21 ) 


where V^ = V^ — V^ ■ According to our assumptions the 
impedances ZjA) are purely imaginary, i.e. Reil/Zj) = 
0. That is why the terms oc V~{t')V~{t") do not appear 
in iSem- In contrast, such terms present in the action of 
the resistors (ITSll even if one puts Ai = A 2 = 0. These 
terms are related to dissipation in the resistors and de¬ 
scribe the current noise associated with it. 
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At long observation time, t ^ 1/?}, l/wc, the full ac¬ 
tion © acquires the form 

.s' = IV Ck), ( 22 ) 

2 ILJn 


where a;„ = 2TTn/t are discrete frequencies, V'^{uJn) = 
{V^{uJn),V^{uJn),V^{uJn),V^{uJn)) is the vector of 
Fourier transformed voltages, and 


MxU^n) = 


/ 

2r,++iui„) 

1,1 'nt 

0 

Ri 

Zq Ri 


1 1 

2 r;, {uj„) 

1 


Z\ Zq R\ 

0 

Ri 

1 

Zn 


^ 0 * 

R 2 

V 

1 

^0 

0 

1 1 vt 

Z 2 Zq 
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(^n ) 
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I)+ i-Uln) 


0 

2'?r 






(23) 


fl2 


The Gaussian path integral over Vn is evaluated 
exactly. Utilizing the property M\{uj) = 
in the long time limit we find CGF in the form 
= -t/“|^ln[detM(w)/detMA=o(w)]. Evalu¬ 
ating the determinants, and keeping in mind that Z* = 
— Zj for reactive elements, we find 

poo j 

F{t,X) = -t J ^ In [l - r(a;){ni(a;)[l-b 712 ( 0 ;)] 

X (e*^*^ — 1 ) -b [1 -b r7i(o;)]772(o;) (e ~ l) }] ■ (^4) 


Here A = Ai — A2, 


I \UJ I I 

RiR2\Gi + G2 + ZoGiG2\ 

is the effective transmission probability, and Gj = 1 /Rj -b 

l/Z,{u). 

Equation (|24ll is the main result of our paper. It 
is the CGF of photons which are scattered by a bar¬ 
rier with the transparency t{uj) and carry the energy 
uj each. It is consistent with standard results of quan¬ 
tum opticsi^ and closely resembles the CGF of scattered 
electrons^, which are fermions. In the context of pho¬ 
ton scattering by a cavity similar expression has been 
derived by Beenakker— , and in the context of phonon 
heat conductance — by Saito and Dhar— . If both ni(o;) 
and 712(w) have the equilibrium Bose form, CGF (l24l) 
acquires the property F{X) = F{—X + 
which translates into the fluctuation theorem P{Q) = 
P(—Q) exp[(5(T]“^ —T 2 "^)]. We remind that the Eq. (1^ 
has been derived assuming Gaussian fluctuations of cur¬ 
rents and voltages in the electric circuit. That implies, 
in particular, that the resistors Ri and R 2 are linear ele¬ 
ments, which do not exhibit Coulomb blockade or other 
types of non-linearities. Besides that we have assumed 
that the real parts of the impedances Zj{uj) are equal to 
zero and they correspond to purely reactive elements like 
inductors, capacitors or their arbitrary combinations. 
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FIG. 2. Distribution of energy transmitted between the re¬ 
sistors during time t for different transmission probabilities 
(a) t{ui) =const, Ti = 300 mK, T 2 = 100 mK, the 
observation time is f = 10 ns. (b) t{uj) has the Lorentzian 
shape, Ti = 300 mK, T 2 = 100 mK, f = 1 ms, CGF is given 
by Eq. (1271) . Discrete number of transferred photons n is 
shown on the horizontal axis. 


III. RESULTS AND DISCUSSION 


Let us now consider some limiting cases. First we as¬ 
sume that the transmission probability, r, is constant 
and the photon distribution functions have equilibrium 
Bose form. In this case the heat current acquires fa¬ 
miliar form Jq = —{i/t)dF/dX\x=o = ttt (Pf — T 2 ) /12. 
The simplest example of such a system is given by two 
directly connected resistors (Fig. le), in which case 
r = 4i?ii?2/(Pi + ^ 2 )^- In Fig. [5^ we show the distri¬ 
bution P{t, Q) for three different values of r. The distri¬ 
bution becomes Gaussian at sufficiently long observation 
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time such that Jqt ^ Ti. The low frequency noise of the 
heat current is given by the expression 


t d\^ 


A=0 


(^r{l-T) + ^T^]{T^+Ti) 


. r 

^ 27r (e-/Ti _ 1) (e-/T2 _ 1) ■ 


(26) 


Another interesting limit is transmission within a nar¬ 
row Lorentzian with r(w) = ri„axr^/[(w — wq)^ -I- T^] and 
r <C Wo, Ti, T 2 . In this case 


F = -n (x/l - w/(wo) - 1 ) , (27) 

where /(wq) = ni(wo)[l n 2 (wo)] - l) -h [1 

ni(wo)]n 2 (wo) (6“*^“^“ — l) . Since F{X) becomes a peri¬ 
odic function of A in this approximation, we get P{t, Q) — 

Y.^Pn5{Q-nuJo) withp„ = ^ dAbe¬ 

ing the probability to transmit n photons with one fre¬ 
quency Wo. The distributions p„ for three different values 
of Tmax are shown in Fig. [21 d. Due to the suppression of 
the average heat current between the resistors the distri¬ 
butions Pn significantly deviate from the Gaussian form 
even though the observation time is long, t = \ ms. It is 
obvious from Eq. (1221) that the distribution becomes 
Poissonian in the limit Ti ^ T 2 and Tmax 1. At higher 
transparencies it deviates from the Poissonian form sim¬ 
ilarly to what has been predicted in Ref. [U, where the 
statistics of photons emitted by a coherent conductor has 
been studied and rectangular shape of the transmission 
line has been assumed. The average heat current and the 
noise corresponding to CGF (l27)) are (here rij = nj{oJo)) 


•^Q =-^rwo[ni - 712], Sq = ^^Fwo(ni[l-hn 2 ] 

-I- [1 -I- ni]n2 -I- "^max [ni-n2]V2). (28) 


Next we assume that leads are attached to the resistor 
1 and bias current / is applied to it (see Fig. [3^). The 
electron distribution function inside it acquires a non¬ 
equilibrium double step formic, f{E,x) = {x/Li)f{E — 
eV) + {l—x/Li)f{E), where V = IRi is the voltage drop. 
We also assume that the temperatures of the resistor 2 
and of the outer leads are much lower than eV. In this 
case one can put 77.2 (w) = 0 and from the Eq. dni) we 
find ni (w) = {eV — uj)9{eV — uj) /6w. Thus the GGF (l24l) 
takes the form 


F = -t 


r-^in 

T 27r 


1 _ I Awi _ 1) 

6 w '' ’ 


■ (29) 



FIG. 3. (a) Bias current I is applied to the resistor 1 in order 

to drive it out of equilibrium. Two capacitors C, which shield 
the detector resistor 2 at low frequencies, are big enough to 
become fully transparent at frequencies w max{ri, T 2 , eV}, 
where V = IR\. In this case the barrier transmission r may 
be approximately treated as frequency independent constant, 
(b) Distribution of transmitted energy during the observation 
time t = 100/eF for three different values of r. Q and P{Q) 
are scaled with the characteristic photon energy eV. 


derived. From CGF 
model 

T{eVY 

2477 


Jq = 


we find these parameters in our 
T{eVf 




727r 


(30) 
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In the classical limit ^ Wc CGF dMl) reduces tc 
F = -t [ p\n[l- r(w) (7AAri2 - X^TiT 2 )] (31) 

Jo ZTT 

where AT 12 = T 1 —T 2 . It is interesting to compare this re¬ 
sult with the experimen t ^^ 4 ^ . In that experiment capac¬ 
itors have been used, which implies Zj{uj) = l/(—fwCj) 
(see Fig. |4^). Accordingly, t(w) (l25l) takes the form 

2/3(wto)^ 


t(w) = 


1 -I- 2(a - l)(wto)^ + (wto)^ ’ 


(32) 


with to — \/ RiR 2 {CiC 2 + CqCi C 0 C 2 ), 

a = 1 + [Rl(Co + Cif + Rl(Co + C 2 )" + 2RiR2C^]/2tl, 

and /3 = 2C'o/(C'i(72 -t-CoCi -I-C 0 C 2 ). For this model one 
can exactly evaluate CGF (ISTI) . 

^ t ( ^(7AAri2 - A2riT2) 

^-To[y2-y2 2 



The corresponding distribution P{Q) is shown in Fig. 
[5b- It is strongly asymmetric with P{Q) = 0 for Q < 0, 
i.e. over long intervals of time, eVt ^ 1, the energy flows 
from the biased resistor to the unbiased one, but never 
in the opposite direction. A somewhat similar system, 
namely a biased resistor coupled to an open transmission 
line, has been earlier considered in Ref. [l^, where the 
average value of the heat current and its noise have been 


and the distribution of the transferred heat P{t, Q) = 
f ^ which reads 


P(t,Q) 


t I 2a 

n\j pTiT2t^ + 2QHl 


e*o 




Q 



( 34 ) 
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FIG. 4. (a) Setup of the experimen t^^'^^ . The circuit pa¬ 

rameters are: R\ — R 2 = 10 Mil, Go = 100 pF, Gi = 680 pF, 
C 2 = 420 pF. The parameters defined in the text take the val¬ 
ues a = 2.134, /3 = 0.0506, and to = 6.29 ms. (b) Distribution 
of energy transmitted during the time t = 0.1 sec and for resis¬ 
tor temperatures Ti = 296 K, T 2 = 88 K. Circles - experimen¬ 
tal pointaiiii^; blue line - Eq. IMI; red dashed line - Gaussian 
approximation P{t,Q) — exp[—(Q — Jqt)^/2SQt]/^2nSQt, 
where Jq and Sq are defined by Eqs. (I35|. 


Here Ki{x) is the modified Bessel function of the second 
kind, and a = a/2 + f3{Ti — T 2 y/ 8 T 1 T 2 . One should bear 
in mind that the expression (1341) is valid in the long time 
limit t ^ to- The average heat current from the resistor 
1 to the resistor 2 and the corresponding noise in this 
model have the form 

m - T2) „ pt,T2 , PHTi-T2)\ ^^, 

^ 4y2a3/2to ^ ’ 

We compare the distribution (1341) with the experimen¬ 
tal oneiiii^ in Fig. 0 )d. The agreement between the two 
is quite good. In particular, one can see the deviations 
from Gaussian form at the tails of the distribution. The 
subtle point of the measurement o^^d^ was the difference 
between the heat Qi, i.e. the change of the energy of 
the resistor 1, and the work ITi, which also includes the 
change of the electrostatic energy of the capacitor Ci. 
We have verified that in the long time limit both Qi 
and Wi should have the same distribution (IMl) . On the 
qualitative level this can be understood from the relation 


Wi = Qi + Ci[Vi{t) — Vi{G)]/ 2 . Indeed, the average 
value of the last term, i.e. of the change in the energy 
of the capacitor Ci during the observation time t, equals 
to zero because {Vi{t)) is finite and does not grow in 
time. Since both Qi and Wi grow in time linearly, one 
can put Wi ~ Qi at sufficiently long t even without av¬ 
eraging. Experimentally, however, the work distribution 
has approached the long time limit form faster than the 
heat distribution. That is why in Fig. |41 d we plot the 
experimental work distribution P{Wi). Further analy¬ 
sis is required in order to understand the origin of this 
behavior. 

We propose the distribution of heat in the low temper¬ 
ature quantum regime to be measured in the setup sim¬ 
ilar to the one used in the experiments [ii^ . Namely, 
one would monitor the temperature of the detector resis¬ 
tor 2 in real time with the time resolution of the order 
of to ~ 12 ‘h/TTTkBTi^ that is the time interval during 
which an average energy ksTi is transferred from the 
resistor 1 to the resistor 2. Assuming Ti = 100 mK 
and r = 3 X 10“^ one hnds to ~ 1 MS, which is within 
the reach of current technology^^. The expected magni¬ 
tude of temperature fluctuations in the second resistor 
caused by fluctuations of heat flow may be estimated 
as <5T| Ri 'iTT 2 t/TT^Tiksv^^^, where t is the observation 
time. For a resistor with the volume ^2 = 0.001 p,m^ 
made of copper (density of states v r: 10 ^® ^m“^) 

and for T 2 = 50 mK and t = lOOto one finds 6 T 2 ~ 15 
mK, which is measurable with currently available ther¬ 
mometers based on normal metal - superconductor tun¬ 
nel junction o^^i^^ . One can further optimize the system 
by, for example, designing the coupling circuit with nar¬ 
row line transmission spectrum, or by using other types 
of temperature sensors like, e.g., recently proposed sensor 
based on an SNS Josephson junctio n^^i^^ . 

In summary, we have developed a theory of full count¬ 
ing statistics of heat exchange between two metallic resis¬ 
tors, which is valid both at high and at low temperatures, 
where the classical formula for the noise © can no longer 
be used. Fluctuations of the heat current in this system 
can be interpreted as scattering of photons by an effective 
potential barrier. In high temperature limit our results 
are in good agreement with recent experimen t ^ ^ 1 ^^ . We 
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